Bose-Einstein condensation in a two-dimensional trap 
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The theory of Bose-Einstein condensation in a two-dimensional(2D) harmonic trap is developed 
from 2D Gross-Pitaevskii equation. The 2D interaction strength is obtained from a 2D collision 
theory. We show the realization of 2D condensation of trapped Bose atoms directly by obtaining 
the stable solutions for the condensate wave function from the 2D Gross-Pitaevskii equation. We 
calculate the ground-state energy of the 2D system, and also the wave function of the 2D vortex 
state. In particular, the 2D energy state becomes less stable than the 3D case with the number 
of trapped atoms. The results of the 2D Bose condensation are also compared with those of the 
well-known 3D case. 
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I. INTRODUCTION 

The recent experimental realization of Bose-Einstein 
Condensation(BEC) of alkali atoms in magnetic traps 
|^|-^ has generated much interest and activity in theo- 
retical and experimental physics. An interesting, but less 
addressed, question is whether such a phase transition 
due to global coherence exists in two-dimensional(2D) 
space. Even though it has been known that 2D BEC of 
a uniform(untrapped) Bose gas cannot occur in 2D at a 
finite temperature since thermal fluctuations destabilize 
the condensate , it has been suggested that spatially- 
varying potentials which break the uniform distribution 
may create BEC in 2D inhomogeneous systems P-[ll| . In 
the presence of trapping potential, the effect of thermal 
fluctuations are strongly quenched due to the different 
behavior exhibited by the density of states. 

At zero temperature, the 2D condensation can be de- 
scribed by the 2D Gross-Pitaevskii equation(GPE), a 
mean-field approximation for the macroscopic wave func- 
tion of weakly interacting bosons. The GPE is obtained 
from mean-field many-body quantum-statistical theory 
and has been proven to explain the condensate state of 
dilute Bose system satisfactorily. Although no direct ex- 
perimental observation of 2D BEC has been reported yet, 
recent theoretical works suggest a possibility of 2D BEC 
in a trapped condition. However, most of the theoretical 
approaches for 2D BEC have been rather indirect and 
sometimes even faced with several difficulties. 

Tempere and Devreese studied harmonically- 
interacting bosons in 2D. They calculated the critical 
temperature from a grand-partition function to show the 
occurrence of 2D BEC, but without any knowledge of 
the 2D condensate wave function for the density pro- 
file. Jackson et al. ^ studied the 2D vortex state by 
direct substitution of the 3D interaction strength g^u = 
4:T:h^a/m (a is the s-wave scattering length in 3D and m 
is the atomic mass) into the 2D GPE, which resulted in a 
dimensional inconsistency. Bayindir et al. M solved the 



2D GPE using the two-fluid model a density estimation 
and obtained the temperature dependence of the internal 
energy and condensate fraction, where the 2D interaction 
strength was still treated just as a free variable. Haugset 
et al. 10 calculated the density proflle and ground state 
energy for a flnite number of 2D Bose particles by diag- 
onalizing the Hamiltonian numerically, where a modified 
2D interaction strength was used to solve the dimensional 
inconsistency. Analytical approach of interacting Bosons 
in 2D trap by 2D nonlinear GPE was attempted by Gon- 
zalez et. al. |TT| , and they studied the ground-state en- 
ergy density by applying the Haugset et aZ.'s modified 
form of the 2D interaction strength to the 2D GPE. 

Previous 2D researches are based on the size- 
independent concepts: if the atomic motion in any one 
axis is frozen completely in the ground state, the sys- 
tem is considered as two-dimension. Moreover, most 
of the studies considered the 2D trap as an extremely 
anisotropic 3D trap and still far from the pure 2D system. 
On the other hand, we will focus on the 2D system which 
is confined by a 2D harmonic trap in (a;, 2/)-direction and 
by ideal rigid walls in ^-direction ||l2| . In this system, the 
2D interaction strength has a logarithmic form which will 
be derived using a 2D scattering theory in the next sec- 
tion. Quasi-2D atomic systems, where atoms follow the 
2D kinematics with the 3D interactions [|l3|, have been 
recently realized for spin-polarized atomic hydrogen ad- 
sorbed on the surface of liquid helium 0,11^ and quasi- 
condensation has been indirectly demonstrated Al- 
though no pure 2D BEC has been studied experimentally 
yet, the 2D GPE is quite valid even in quasi-2D at low 
temperatures . Moreover it will provide us a guide to 
a better understanding of the 2D statistical properties of 
the quasi-2D system. 

In this paper, we will solve the 2D GPE numerically 
and show the existence of the 2D BEC directly by ob- 
taining the stable condensate wave functions and the 2D 
ground state energy density of the trapped Bose atoms. 
Possible vortex states in 2D are also discussed. The re- 
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suits of the 2D Bose condensation are then compared 
with those of the weU-known 3D cases. 



II. TWO-DIMENSIONAL SCATTERING THEORY 

The relation between interaction strength, g, and the 
s-wave scattering length, a, is well-known in 3D as g — 
47r ?i^a/m. However, the 3D result is not applicable to 
a 2D GPE: substitution of the 3D relation into the 2D 
GPE results in dimensional inconsistency. As a first step 
to obtain the correct 2D GPE, the interaction strength 
in 2D is derived from the following 2D scattering the- 
ory [p7|-pO[. We begin with the 2D time- independent 
Schrodinger equation with the interaction U{p), given by 



u{p)Mp). 
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where fc^ = 2/i„i_E/?i . is the reduced mass of two 
particles (m/2). One can find the general solution of Eq. 
(0) with the help of the 2D Green's function 
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The general solution of the 2D Green's function is the 
Hankel function |^,^ . Therefore the 2D wave function 
can be expressed as 
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Here Hq is the first kind Hankel function of order zero 
defined by Hq(x) = Jq{x) + iiVo(x), where Jq and 
Nq are Bessel and Neumann function of order zero. 
For large p, it has the asymptotic behavior Hi{x) ~ 
^^/2y^gi(a:-i7r/2-7r/4)^ where I is an integer. The Born 
approximation is applied to the last step, and k' = kp. 

In 2D, for large p, the scattered wave function should 
have the form 



Fk- 
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where Fk is the first-order Born scattering amplitude in 
2D which has a dimension of {length)-^ ^ . This asymp- 
totic formula should be a solution of the time indepen- 
dent Schrodinger equation (0), and satisfy p — s- cx) limit 
of the Green's function solution. The scattering cross- 
section in 2D now has the dimension of length given by 
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where jinc and jsc are incident and scattered flux densi- 
ties. We note that the total scattering cross-section for 



small k should be 27r|Ffcp in 2D instead of 47r|/3£ip in 
3D. If one assumes the delta- function type interaction 
U{p) = {mg/h^)6'^{p), the 2D scattering amplitude Fk is 
obtained in the complex form from Eqs. (H) and M) as 
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The next step is to find the relation between the scat- 
tering amplitude and the scattering length in 2D. The 
well-known relation fsjj — —a, where f^D is the 3D scat- 
tering amplitude and a is the s-wave scattering length 
in 3D, is not directly applicable in 2D. In general, the 
s-wave scattering length in 2D, b, is different from the 
one in 3D, a, and the s-wave scattering length in 2D is 
not known yet from experiment. 

The relation is obtained from partial wave analysis of 
collision theory. The incident wave can be written as an 
expansion of plane wave in cylindrical coordinate l22] 
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where I is an integer. Away from the range of potential, 
the scattered wave has a phase shift Si , and the scattering 
matrix is obtained as Si — e^'*' since = 1 for elastic 
scattering. Therefore, for large p, the solution of Eq. (|l|) 
is written as 

-. oo 

Mp) = e^"^-" + 2 E *'(^' - miikp)e'"' 

1 — — 00 

^ oo 

:^e'^-P + ^^= y i^{e^'^' -l)e'^''P-^-^^e'^'^. (8) 
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Scattering length is defined as the distance where the 
two-body wave function vanishes for zero energy. The 
phase shift due to potential scattering is given as a func- 
tion of the scattering length, and it is well-known as 
'5o,3Li = —ka in 3D, whereas 
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in 2D [Q. Note that it is effective only in the low en- 
ergy scattering limit, fc6 << 1. It can be easily checked 
as follows. The scattering length is just the intercept of 
the radial wave function satisfying the boundary condi- 
tion = 0. Therefore, for / = 0, the wave function 
at large distance and small k is expressed as 

'''' --y^h 
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The necessary condition for the vahdity of the Born ap- 
proximation is that the phase shift So be very small for 
small k, which can be easily confirmed from Eq. (H). 
Note that unlike the 3D case where a can be negative, 
we do not consider the negative scattering length here 
since the centrifugal potential of the lowest partial wave 
is negative in 2D so that the extrapolated local wave 
function cuts the radial axis always above the origin . 

Comparing Eq. (^) with Eq. (Q), the scattering am- 
plitude is written as a series expansion, and the I = 
state contributes to the 2D system. Therefore, the 2D 
scattering amplitude becomes 
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Finally, we obtain the 2D interaction g from Eqs. 
and ([n] ) as 
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In general the 2D interaction strength is given as EtI 
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Here ^ is a dimensionless atomic parameter given by 
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where ^ is positive since fc6 <C 1. Note that the logarith- 
mic dependence of the interaction term suggests similar 
condensate characteristics in 2D for most bosonic alkali 
atoms with positive scattering length. 

In addition to the 2D scattering length, b, we do not 
have any reliable value of the wave-vector, k, in 2D. An 
approximation that the wave-vector k be the inverse of 
the largest distance available in the perpendicular direc- 
tion may not be correct. However, the experimental re- 
sults provide the value of the product kb as described in 
the next section. Note that in an extremely anisotropic 
3D system, the dimensionless atomic parameter is given 
as = (l/V27r)a/az, where Qz = ^Jhjmwl, and more- 
over < C ill general. 



III. THE CONDENSATE STATE AND ITS 
ENERGY STATE 

Now the 2D condensate wave function of trapped di- 
lute Bose atoms of mass m can be obtained from the 2D 
GPE 



~vl^v,M^N9r{p) 
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V'(p) = /xV'(p) , (15) 



where / (Pp'4P{p) = 1, and A'' is the number of conden- 
sate particles. Here, we assume a 2D isotropic, harmonic 
trap Vext{p) = \rnijj^p^ where uj is the trap frequency. 



Then we can simplify Eq. ( p_5D for numerical calcula- 
tion by introducing dimensionless variables (p — > auoP, 
p — > hujp, and -0 — > fl/^JV') 

[-V2 + p2 _ 2^ + 8^e^^'(p)] ^(P) = 0, (16) 



where auo — ^/Ti/mu) is the 2D harmonic oscillator 
length, and p is the 2D chemical potential which is ob- 
tained from the normalization condition. Note that the 
product of incident wave-vector and scattering length, 
kb, is the only atomic parameter that contributes to the 
condensate states. 

In the numerical calculation of the GPE, we need the 
value of the atomic parameter ^ in Eq. (nSh which is a 
function of k and b. Although they are not known sep- 
arately, we are able to deduce their product kb from ex- 
perimental data, which appears together in the 2D GPE. 
It is certain that different Bose atoms have different 2D 
scattering lengths, but the logarithmic dependence makes 
the difference less sensitive. A recent experiment of the 
hydrogen on helium surface by Safonov et al. ||l5[] has re- 
ported ^ = 1/7 which indicates kb = Zx . Although 
their system actually satisfies a quasi-2D condition, one 
can still take this value in the 2D GPE as an effective 
interaction potential. Since it is well known that the 2D 
GPE is valid even in quasi-2D at low temperatures , 
and moreover it hints the 2D atomic characteristics in 
the quasi-2D system, it is useful to quantify the criterion 
for the quasi-2D, or the effective thickness of the 3D sys- 
tem to exhibit the 2D statistical properties. The criteria 
of the effective thickness of the trapped 2D Bose system 
may by obtained |l^ . 

The next procedure to solve Eq. ( p^ ) is almost similar 
to those of 3D We have plotted the 2D conden- 

sate wave functions versus p for several values of atom 
number N in Fig. 1. It corresponds to the z = cut 
of the anisotropic contour plot of the ground states in 
3D. The spatial distribution of the condensate in 2D is 
much broader than that in 3D, and the condensate wave 
functions approach the parabolic limit more rapidly with 
the increase of the number of atoms. In other words, 
the effect of atomic interaction potential becomes more 
prominent in 2D. 

In the non-interacting case, the solution is still Gaus- 
sian with Tpip) = 7r~^/^e~'' In the strongly repulsive 
limit or Thomas- Fermi limit, it has the parabolic solu- 
tion of V'^(p) = (2m — p'^)/8t^^N. The overall shape of 
the condensate wave functions are similar with that of 
3D, but it approaches the parabolic limit more quickly 
with the number of atoms Therefore, the peak of 

the density profile decreases much faster in 2D. The 2D 
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healing length that balances between the quantum pres- 
sure and the interaction energy of the condensate is also 
different from that of 3D. Refer to TABLE 1 for detailed 
comparison. 

The ground-state energy for 2D condensate bosons can 
be calculated in a similar way. With the dimcnsionless 
variables defined before, we obtain the dimcnsionless en- 
ergy density as 



Sp. (17) 



Using a Gaussian trial function, we easily find the 
ground-state energy density satisfies 



(18) 



The ground-state energy per particle in Eq. (|18|) is plot- 
ted in Fig. 2, and compared with the well-known 3D re- 
sults of ^^Rb. The 2D system becomes more excited for 
a given N and less stable as N is increased with respect 
to the 3D case. We have summarized the fundamental 
differences between our 2D results and the well-known 
3D ones in TABLE 1. 



IV. THE VORTEX STATE 

Now let us consider the vortex states of the 2D system. 
The hydrodynamic theory connected to superfiuidity is 
understood by vortex. The 2D system can be rotating 
about the center of the 2D trap to give quantized circula- 
tion of atomic motion. The angular momentum quantum 
number n gives the quantum winding of the 2D vortex 
state. A vortex state with winding number k is written 
as V'(p) = ^(p)e''^(''\ where (f>{p) — \/n{p) is the modu- 
lus. The phase S is chosen as k6 where k is an integer. 
An angular momentum quantum number k can be then 
assigned to the quantum winding of the 2D vortex state, 
and one finds the vortex states with a tangential veloc- 
ity V — Kh/mp. As a result of the quantum circulation, 
the angular momentum of the system with respect to the 
p = axis becomes L — Nnh. 

Adding the vortex term of K^/p^, Eq. (16) is directly 
converted into the vortex state 
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The wave function for k = 1 vortex state is plotted in Fig. 
3. The overall shape and N dependence of the vortex- 
state wave function are similar with those of 3D . As 
expected, the vortex state also corresponds to the z — Q 
cut of the anisotropic contour plot of vortex state in 3D. 
We also observe that the 2D vortex has a larger radius 
than the 3D one. 

The critical angular velocity or the energy difference 
between the vortex state of k = 1 and the ground state 
K = in Eq. dll) is obtained analytically as 
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(20) 



For £k=i, the trial function of 



pe- 



type was em- 



ployed. With the increase of the number of particles, the 
vortex-excitation energy in 2D is much smaller than the 
3D one, which indicates that vortices are expected to be 
produced more easily in 2D than 3D. 

Although there are fundamental differences between 
the 2D and the 3D results as summarized in TABLE 
1, it will be interesting to compare with the quasi-2D 
scheme considered as a limiting case of 3D. For compar- 
ison of the 2D trap with an extremely squeezed 3D trap, 
we take the 3D external potential given by Vextip, z) = 
{\ / 2)muP' {p^ ~V X'^ z^) . Here A is the anisotropy parameter 
and much larger than 1 for the quasi-2D system, whereas 
X oo for the 2D trap. With a fixed number of atoms, 
we plot the condensate wave functions versus A, from the 
3D GPE in Fig. 4. We find the 3D wave functions merge 
very slowly with A, but does not reach the pure 2D limit 
at all. 



V. DISCUSSIONS 



The 2D nonlinear Schrodinger equation(GPE) is not a 
simple extension of the 3D case but is connected into a to- 
tally different 2D collision theory which requires different 
approach. We have developed the theory of the 2D GPE 
for trapped neutral Bose atoms in a 2D harmonic trap. 
Applying the quasi-2D experimental value of ^ = 1/7 to 
the effective interaction strength, one can solve the 2D 
GPE numerically without detailed knowledge of k and h 
separately. 

We have obtained the stable solutions of the 2D 
GPE, which predicts possible existence of 2D BEG. The 
ground-state energy of the condensate particle is also cal- 
culated and it is found that the 2D system becomes less 
stable as the number of trapped atoms is increased, com- 
pared with the 3D case. We also have obtained the wave 
functions of the 2D vortex state numerically and its crit- 
ical angular velocity. 

The 2D BEG transition may look similar to the 
Kosterlitz-Thouless(KT) |26) vortex-state transition, but 
the phase transition of 2D BEG does not need any strong 
interaction between atoms. That is the fundamental dif- 
ference between the two transitions. 

The logarithmic dependence of the interaction poten- 
tial on the scattering length makes the 2D system less 
sensitive to the species of the atoms used for condensa- 
tion. That will make every alkali atoms of positive scat- 
tering length show similar condensate wave patterns in 
2D. The comparison between our 2D and well-known 3D 
results are summarized in TABLE 1. 

The quasi-2D scheme can be interpreted as a limit- 
ing case of the 3D one. By varying the trapping field. 
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it is possible to separate the single-particle states in the 
oscillation into well-defined bands. To compare the sit- 
uation of 2D trap with an extremely squeezed 3D trap, 
we take the external potential for 3D as Vextifj^, z) — 
{l/2)mLj^{r\ + A^z^), where A 1. The case of nega- 
tive scattering length in 2D will be discussed later. 
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FIG. 1. The 2D condensate wave function of neutral 
trapped bosonic atoms for kb — 3 x 10"^. 



FIG. 2. The ground-state energy per particle (in unit of 
huj) of Eq. (10) for kb = 3 x 10"^ 



FIG. 3. The 2D wave function of vortex state of neutral 
trapped bosonic atoms for kb — 3 x 10~^ and a quantum 
circulation k = 1. 



FIG. 4. The condensate wave functions in anisotropic 3D 

trap for N = 100. From up to down A = 2, 10, 50, 100, 500 . 

,1/2 



Note that 0(p) = [/ dz\4}{p,z)\^ 



TABLE I. The comparison of the 3D and the 2D BEC. The dimensionless atomic variable ^ is defined by ^ = — l/(21nfc6). 
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Dimension 


D = 3 


D = 2 


Scattering amplitude 


ho = -a 


F2D = 


Interaction strength 
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Healing length 
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Chemical potential 
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